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Abstract 



Based on the assumption that D* 7 (2317) and Z? s j(2460) are the (0 + , 1 + ) chiral partners 
of D s and D*, we evaluate the strong pionic and radiative decays of D* sJ (2317) and Z? s ,/(2460) 
in the Constituent Quark Meson (CQM) model. Our numerical results of the relative ratios 
' of the decay widths are reasonably consistent with data. 

^ PACS numbers: 13.25.Ft, 12.38.Lg, 12.39.Fe, 12.39Hg 

■ 1 Introduction 



The new discoveries of exotic particles Z)*j(2317) and D s j(2460), which possess spin-parity 
structures of + , 1 + respectively ^HHH], attract great interests of both theorists and experi- 
mentalists of high energy physics. Some authors |3] suppose that D*j(2317) and D s j(2460) are 
(0 + , 1 + ) chiral partners of D s and D* i.e. p-wave excited states of D s and D* [H]. Narison used 
the QCD spectral sum rules to calculate the masses of D*j(2317) and D s j(2460) by assuming 
that they are quark-antiquark states and obtained results which are consistent with data within 
a wide error range (Sj. Beveren and Rupp also studied the mass spectra and claimed that 
their results support the cs structures for D* sJ (2317) and D s j (2460). Meanwhile, some other 
authors suggest that D*j(2317) and D s j(2AQ0) can possibly be four-quark states |8l 19} ITOl ITT] . 
Thus, one needs to try various ways to understand the structures of D*j(2317) and Z? s j(2460). 
In general, one can take a reasonable theoretical approach to evaluate related physical quantities 
and then compare the results with data to extract useful information. One can determine the 
structures of D*j(2317) and -D s j(2460) by studying the production rates of the exotic particles, 
and our recent work ^2] is just about the production of D*j(2317) in the decays of ^(4415). 

Recently, several groups have calculated the strong and radiative decay rates of L>*j(2317) 
and D s j(2460) in different theoretical approaches: the Light Cone QCD Sum Rules, Constituent 
Quark model, Vector Meson Dominant (VMD) ansatz, etc. [THl 1131 H3 [TH1 HZJ [TH1 HH1 - For a 
clear comparison, the results by different groups are listed in Tables 2 and 3. The authors of Ref. 
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[§1 I2(jj also calculated the rates based on the assumption that D*j(2317) and D s j(2460) are in 
non-cs structures. Their predictions on the D*j(2317) and D s j(2460) decay rates are obviously 
larger than that obtained by assuming the two-quark structure by orders. Thus studies on 
the strong and radiative decays with other plausible models would be helpful. It can not only 
deepen our understanding about the characters of these particles, but also test the reliability of 
models which are applied to calculate the decays. Because D s j(2632) was only observed by the 
SELEX collaboration [2H, but not by Babar [221, Belle [23] and FOCUS [21, its existence is still 
in dispute, so here we do not refer decays of D s j(2632). 

In this work, we study the strong and radiative decays of Z)*j(2317) and -D s j(2460) in the 
framework of Constituent-Quark-Meson (CQM) model. CQM model was proposed by Polosa 
et al. [25] and has been well developed later based on the works of Ebert et al. [22] (See the 
Ref. |22| for a review). The model is based on an effective Lagrangian which incorporates the 
flavor-spin symmetry for heavy quarks with the chiral symmetry for light quarks. Employing 
the CQM model to study phenomenology of heavy meson physics, reasonable results have been 
achieved [23 [2E]- Therefore, we can believe that the model is applicable to our processes and 
expect to get relatively reliable conclusion. 

The constraint from the phase space of the final states forbids the processes D*j(2317) — > 
D s rj(r]') and D s j(2AQ0) — ► D*rj(rj'), so that the only allowed strong decay modes are D*j(2317) — ► 
D s tt° and L> sJ (2460) -> D*n°. In principle, D sJ (2AQ0) -> L>*j(2317) + vr° which is a 1+ 
+ + 0~ process, is allowed by the phase space. However, it is a p-wave reaction, and the total 
rate is proportional to |p| 2 , where p is the three- momentum of emitted pion in the center-of- 
mass frame of D s j(2460). In this case |p| is very small (about ~ 28 MeV), so that this process 
can only contribute to the total width a negligible fraction, in practice. 

The aforementioned strong decay modes obviously violate isospin conservation. Therefore 
the decay widths of £^(2317) -> D s tt° and £> sJ (2460) -> D*tt° must be highly suppressed. 
Moreover, direct emission of a pion is OZI suppressed |29j . 

Cho et al. suggested a mixing mechanism of rj — tt° where the isospin violation originates 
from the mass splitting of u and d quarks |3(J| . In that scenario, L>*j(2317) and -D s j(2460) 
firstly transit into D s rj and D*r], and then rj transits into 7r° by the mixing. In the intermediate 
process, rj obviously is off-shell. The mixing depends on the mass difference of rj and ir, and the 
effects due to the mixing between rf and ir can be ignored. 

Another sizable mode is the radiative decay. Even though, by general consideration, the 
electromagnetic reaction should be much less important than the strong decay, it does not suffer 
the suppression of isospin violation, therefore one may expect that it has a comparable size 
to the strong processes described above. The relevant decay modes are D*j(2317) — > D* + 7, 
D sJ {2Am) -^D s + 1 and £> sJ (2460) -»• D* + 7 and D sJ (2460) -» £>*j(2317) + 7. 

Currently the Babar and the Belle collaborations have completed precise measurements on 
the ratio of r(£> sJ (2460) -» D s j) to T(D sJ (2A60) -» D*7r°) jSJEUE2]. Now Babar and Belle 
collaborations begin to further measure other decays of D*j(2317) and D s j(2460). We are 
expecting new results of Babar and Belle which can be applied to decisively determine the 
structures of £>*j(2317) and D sJ (2A60). 

This paper is organized as follows, after the introduction, in Sect. 2, we formulate the strong 
and radiative decays of D* sJ (2317) and D s j (2460). In Sect. 3, we present our numerical results 
along with all the input parameters. Finally, Sect. 4 is devoted to discussion and conclusion. 
Some detailed expressions are collected in the Appendix. 
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2 Formulation 



First, for readers' convenience, we present a brief introduction of the Constituent-Quark-Meson 
(CQM) model [221 • The model is relativistic and based on an effective Lagrangian which com- 
bines the heavy quark effective theory (HQET) and the chiral symmetry for light quarks, 

Lcqm = x[l ■ (id + V)] X + XI ■ A^x ~ m q xx + ^-Tr[d^d^+] + h v (iv • d)h v 



[ X (H + S + if^)h v + h.c] + ^-Jr[(H + S)(H - S)] + ±-Tr[f^}, (1) 



where the fifth term is the kinetic term of heavy quarks with if)h v = h v ; H is the super-field 
corresponding to doublet (0~, 1~) of negative parity and has an explicit matrix representation: 

1 + 4 

where P and P*^ are the annihilation operators of pseudoscalar and vector mesons which are 
normalized as 

(0|P|M(0-)) = VMh, and <0|JP*"| M (1")) = yfEfaf* . 
S is the super-fields related to (0 + , 1 + ), 

S = ^K^Ts - Po}. 

%M_ 

X = il = u i d, s) is the light quark field and £ = e fa , and M is the octet pseudoscalar matrix. 
We also have 



and 



Because the spin-parity of D*j(2317) and £> s j(2460) are 0+ and 1+, thus D* J (2317) and 
D s j(2A60) can be embedded into the S-type doublet (0 + , 1 + ) [2S], whereas D s and D* belong 
to the H type doublet (0~, 1~ ). Then we can calculate the strong and radiative decay rates of 
D* sJ (2317) and L> sJ (2460) in the CQM model. 

2.1 The transition amplitude of D* sJ (2317) -> D s + tt° and D sJ (2460) -> D* + tt° 
strong decays in the CQM model. 

As discussed above indirect D* sJ (2317) —> D s + tt° and D sJ (2A60) — > D* + tt° occur via two 
steps. In Fig. we show the Feynman diagrams which depict the strong processes D*j(2317) — ► 
D s + i] — > D s + 7T° and -D s j(2460) — ► D* + rj — > D* + 7r°. According to chiral symmetry, r\ only 
couples with light quark, thus it can only be emitted from the light-quark leg in FigflJ 

The matrix elements of £>* 7 (2317) -> D s + vr° and Z) sJ (2460) -> D* + vr° are written as 

M[D* aJ (2317) - D s + vr°] = <7r°|£ m<!riwg |7/) 2 ^ 2 (??£> s |£cQAf |£> s *j(2317)), (2) 
M [Z) aJ (2460)) - + vr°] = (A^mixinglv) 2 l _ 2 (t]D*\Ccqm\D s j(2A60)). (3) 
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^(2317) 




Figure 1: The Feynman diagrams which depict the decays of D*j(2317) — ► D s + 7r° or 
L> s j(2460) — ► D* + 7T°. The double- line denotes the heavy quark (c— quark) propagator. 

The mixing mechanism is described by the Lagrangian 1 

ml rh u - rh d 

L-mixinq — /— ~ 77 

which originates from the mass term of the low energy Lagrangian for the pseudoscalar octet 

m 



2 fl 

4(m u + m d ) 



where m q is the light quark mass matrix. The matrix elements of {r]D s \7i.cQM\D*j(2317)) and 
(t]D*\Hcqm\D s j(2460)) will be calculated in the CQM model. 

fhi{i = u, d, s) are the current quark masses. It is noted that in the CQM model calculations, 
the quark masses (m q , m c ) which we denote as m g , are constituent quark masses (23] , whereas, 
for the mixing, the concerned masses which we denote as fh q are the current quark masses |30]. 

According to the CQM model EH], couplings of L>*j(2317) and D sJ (2460) with light and 
heavy quarks are expressed as 

and 

the couplings of D s and D* to light and heavy quarks are 

i + j> 

2 

and 



ZsM Da ~f 5 



1 + 



Z S M D *J 2 , 



2 

where e± and e 2 denote the polarization vectors of Z) s j (2460) and D* respectively. M\ and M 2 
are respectively the masses of Z?*j(2317) and D s j(2460). The concrete expressions of Zjj and 

Here we ignore the mixing of r\ and rf , because the mixing angle 6 ~ — 11° is small and does not much affect 
our results. Therefore we simply assume that rj is rjs and the contribution of r( , as discussed in the text, is 
neglected in our calculations. 
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Zs are given in Ref. [25] as 

(A H + m s 

dh(A s ) 



Z^ = (AH + m s )^M + / 3 (A 5 ), (5) 



Z^ 1 = (A s + m s )-^ + I 3 (A s ), (6) 



iN c f 1 /^ 2 dy 

h{a) = -T^i] exp[-y(m 2 s - a 2 )](l + erf (a^/y)), (7) 

1d7t 4 Ji/a 2 y i 



where erf is the error function. 

Now, we can write out the transition matrix elements as 

(vD s \Hcqm\D* sJ (2317)) 
- (-I) ly /Z S M 1 Z S M D ^-—J (2?r)4 {l 2_ mm + q) 2_ m 2 ]{v . l + As) . 

(8) 

and 

(^|WcoJifP.j(2460)) 

= {-iy^Z S M 2 Z H M D ^-—J ^ (i2_ m 2 )[(Z + ,)2_ m 2 ](l; . i + As) > 

(9) 

where N c = 3. 

Omitting technical details in the text for saving space, we finally obtain 

(2\ A 



(vD s \H C qm\D* s j(2317)) = J Z s M 1 Z H M Ds (- -)—, (10) 

( V D* s \H CQ M\D s j(2m)) = y/z s M 2 Z H M Dt ( - vf)^^' ( n ) 

with 

A = 4(m 2 s m v uj - m 2 m s )0 - m\{O x + uj0 2 ) + m v uj{20 z -0 4 -0 5 - 2C 6 ), (12) 

where uj = v ■ v' = and the concrete expressions of O and Oi are listed in the appendix. 

The decay widths of D*j(2317) D s + tt° and D sJ (2460) -» L>* + vr° read as 

T[ j d: j (2317) - D s + 7T°] 

'W, (13) 



Z s Z H M Da |p'| / rh u -m d 



10247rM!/2 Vm s - (m u + m d )/2 
T[D S j (2460) D* + vr°] 

ZsZ H M Da |p'| / fhu-fhd 



3072^Af 2 /2 W s - (m u + m d )/2 



2 + Mf M*. 



I^l 2 , (14) 



where the relations = 2mBo and = |(m + 2m s )BQ with m = (fh u + rhd)/2 are employed 
to derive isospin suppression factor (m u — m<i)/[m s — (m u + rhd)/2] [35] . These relations are valid 
at the leading order of the chiral theory, but the principle does not really apply for estimating 
the mass of rf (or r?o) due to an extra contribution from the axial anomaly |33j . 
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2.2 The transition amplitude of _D*j(2317) and D s j(2460) radiative decays in 
the CQM model. 

the Feynman diagrams of D*j(2317) D* + 7, A,j(2460) -> D fl + 7, L> sJ (2460) -> D* + 7 and 
-D s j(2460) — ► D*j(2317) + 7 are presented in Fig. [21 We will evaluate these radiative decays in 



7^9 




£* 7 (2317) 



317) A ! ; + 
2 +P 



(a) 



£>! J (2317) 




Figure 2: The Feynman diagrams depicts the radiative decays of -D*j (2317) and Z^j (2460). 
the CQM model. 

(a) £>;j(2317) ->£>;+ 7. 

In contrast with the case in Fig^ the heavy quark also couples to 7. Thus the transition 
matrix element of D*j(2317) — > D* + 7 is written as 

M[D* aJ (23l7)^D* a +<y] 



(2tt) 4 (I 2 - m 2 s )[(l + q) 2 - m 2 ](v ■ I + A s ) 



2e d 4 Z Tr[(i + m^ 2 ^ 7 i^] 

3i (2vr) 4 (Z 2 - m 2 )(v ■ I + A H )(v ■ I + A s ) 



Z s M 1 Z H M D *(^-){(e~ f -e 2 ) 



+{e 1 -p'){e 2 



m s TZ + Ki + 2TZ 5 



M D * 



(15) 



(b) £> sJ (2460) ^D s + 1 . 

The transition matrix element of -D s j (2460) — > D s + 7 is 



M[Z> aJ (2460)-^Z) a + 7] 



(-l^yZsMaZ^Mz^ ( — 



r«» d 4 Z Tr[(i + m s )^(l + ^ + m s +ie) 75 ^75A] 



(2vr) 4 (Z 2 - m 2 ) [(Z + g )2 - m 2 ] (y ■ Z + A s ) 



2e ri 4 Z Tr[(/ + m s )7 5 ^ 7 ^75A] 1 

3 J (2vr) 4 (I 2 - m 2 )(v ■ I + A H )(u ■ I + A s ) J 



ZsMiZhMd, { ( e 7 ' e i) 



m^-^ 3 -2^-^ ( ^ + ^ + 2 ^ ) 



+(e 7 -p')(ei 



im s n + TZi + 2K 5 



(16) 
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(c) D sJ (2AQ0) D* + 7. 

The transition matrix element of D s j (2460) — ► D* s + 7 is 



M[D aJ (24GQ) -^D* + j} 



{-l)i^Z s M 2 Z H M Dl N c { (-^j J" 
2e r 9 d A l Tr[(/ + m s )^ 2 ^ 7 ^75^i] 



(2vr) 4 (I 2 - m 2 )[(l + q) 2 - m 2 ](v ■ I + A s ) 



3 J (2vr) 4 (P - m 2 )(v ■ I + A H )(v ■ I + A s 
Z S M 2 Z H M D * n±\ EaPpaq ^e^{el 



l<.m s — /Ci H — h 



6 



+ (2^ 5 + 2^ 1 )(p , - (? ) 



M 2 



M Dt M D . 

(17) 



(d) £> sJ (2460) 2^(2317) + 7. 

The transition matrix element of -D s j(2460) — > D*j(2317) + 7 is 



M[D aJ (2460) ->Z£j(2317) +7] 

= (-l),V^1^2iV C {(^) / " "LV-M^-r^....^ 



(2tt) 4 (Z 2 - m 2 )[(Z + g) 2 - m?](c • I + A s ) 



+ 6j (2tt)* (I 2 - m 2 ){v ■ I + A s )(v ■ I + A s 
= VZ S M, Z S M 2 ^H^s^qy^ . (18) 

The concrete expressions of 1Z and IZi are listed in the appendix. 



3 Numerical results 

With the formulation we derived in last section, we can numerically evaluate the corresponding 
decay rates. Besides, we need several input parameters for the numerical computations. They 
include: f n = 132 MeV, m s = 0.5 GeV, A = 1.25 GeV, the infrared cutoff /i = 0.51 GeV 
and A s - A fl = 335 ± 35 MeV m v = 547.45 MeV, Mi = 2317 MeV, M 2 = 2460 MeV, 
Md s = 1968 MeV and Mp* = 2112 MeV [HI]- The suppression parameter was estimated in 
ref.[3S] as 

rh u -rh d 1 
rh s - (m u + rhd)/2 43.7 

(a) We present the decay widths of D* sJ (2317) -> D s +vr° and £> sJ (2460) -> D* + vr° in Table 

□ 

For a comparison, we also list the values of the decay widths of -D*j(2317) — > D s + ir° and 
D s j(2460) — > D* + 7T° which are calculated by other groups, in Tabled 
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Atf(GeV) 


A 5 (GeV) 


Z H (GeV)- 1 


Z s (GeV)- 1 


Ti(keV) 


r 2 (keV) 


0.5 


0.86 


3.99 


2.02 


3.68 


1.86 


0.6 


0.91 


2.69 


1.47 


5.36 


2.72 


0.7 


0.97 


1.74 


0.98 


8.71 


4.42 



Table 1: The values of As and Ah are taken from Ref. According eqs. © and ©, one gets 
the values of Zs and Zh- T\ and r 2 denote respectively the decay widths of Z?*j(2317) — ► D s +tt° 
and £> sJ (2460) -► D* + vr°. 





CQM model 


na 




.16 


El 


m 


Ti(keV) 


3.68 ~ 8.71 


6±2 


34 ~ 44 7 ± 1 


21.5 


~ 10 


16 


r 2 (keV) 


1.86 ~ 4.42 




35~51 7±1 


21.5 


~ 10 


32 



Table 2: In this table, we list our calculation of D*j(2317) -> D s + tt° and £> s j(2460) -> L>* + vr° 
and that obtained by other groups. Where T\ and r 2 are respectively denote the decay widths 
of D* sJ (2317) -> D s + tt° and D sJ (2460) -> D* + vr°. 

Our values depend on the model parameters, but qualitatively, the order of magnitude is 
unchanged when the parameters vary within a reasonable region. All the values in tableware 
somehow consistent with each other for order of magnitude, even though there is obvious differ- 
ence in numbers. 

(b) With the same parameters, we obtain the radiative decay rates of L>*j(2317) and D s j (2460). 
The results are listed in Table. |3] 





CQM 


model 








El 


F(D* sJ (2317)^D s + 1 )(keV) 




1.1 


4~ 6 


0.85 


1.74 


1.9 


r(D sJ (2460) ->£> s + 7 )(keV) 


0.6 < 


- 2.9 


19 ~ 29 


3.3 


5.08 


6.2 


T(D sJ (2460) -Z£ + 7 )(keV) 


0.54 


~ 1.4 


0.6 ~ 1.1 


1.5 


4.66 


5.5 


T(D sJ (2A60) -> £^(2317) + 7 ) (keV) 


0.13 ~ 0.22 


0.5 ~ 0.8 




2.74 


0.012 



Table 3: The rates of D*j(2317) -> D s + 7, L> sJ (2460) -> D s + 7, D sJ (2460) ->•£>* + 7 and 
D s j(2460) — ► D*j(2317) +7. Meanwhile we also list the results obtained by other group. 

For convenience, we define the relevant ratios as 

R x = T(D* sJ (2317) -^D s+1 ): T(D* S ,,(2317) - D s + tt ), 

ii 2 = r( J D sJ (2460) -> D s + 7) : r(D sJ (2460) L>* + tt ), 
i? 3 = r(D sJ (2460) - L>* + 7) : r( J D sJ (2460) -> D* + vr°), 
^4 = T(D sJ (2m) - ££,(2317) + 7) : r( J D sJ (2460) - D s * + tt ). 
Then we calculate the R^s and tabulate the results below. 
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Belle 


Babar 


CLEO 2 


CQM model 


Ri 


< 0.18 31 




< 0.059 


0.12 ~ 0.30 


R2 


0.55 ±0.13 ±0.08 |SI] 


0.375 ± 0.054 ± 0.057 "22] 


< 0.49 


0.32 ~ 0.66 


R3 


< 0.31 31 




< 0.16 


0.29 ~ 0.32 


R4 




< 0.23 32 


< 0.58 


0.05 ~ 0.07 



Table 4: The ratios of radiative decay widths to strong decay widths for L>*j(2317) and 
.D s j(2460). The first three columns are the experimental data and the fourth column is our 
result calculated in the CQM model. 

4 Conclusion and Discussion 

In this work, based on the assumption that D*j(2317) and -D s j(2460) are chiral parters of D s 
and D*, we calculate the rates of £>*j(2317) -> D s ir° and £> sJ (2460) -> D*tt° in the Constituent- 
Quark- Meson (CQM) model and take into account the r\ — 7r° mixing mechanism |30j . We also 
estimate the rates of L>* 7 (2317) -> D* + 7, D sJ (2460) -» D s + 7, D sJ (2460) -> D* + 7 and 
£> sJ (2460) -» 1^(2317) + 7 in the same model. 

Comparing our results about the strong decay rates with that obtained by other groups, we 
find that our results are reasonably consistent with the values listed in table El 

For the radiative decay, our results generally coincide with that obtained by other groups 
and especially these results of the QCD sum rules. 

The ratio R2 has been measured with relatively high precision [3*1 1311 "3*2*j . however for R±, R3 
and R4, there only are upper limits given by Babar, Belle, and CLEO collaborations [2*1 131*1 13*2*1 . It 
seems that our results on the ratios well coincide with the experimental values. This consistency 
somehow implies that the assumption of D*j(2317), D s j(2460 being p-wave chiral partners of 
D s , D* does not contradict to the data with the present experimental accuracy. 

The experimental upper bounds on the total widths are r(D*j(2317)) < 4.6 MeV and 
r(L> sJ (2460)) < 5.5 MeV [23]. Obviously, the overwhelming decay modes of L>* 7 (2317) and 
D s j (2460) are their strong and radiative decays, therefore we can roughly take sums of these 
widths as the total widths of D*j(2317) and D s j(2460). However, our numerical results as well 
as that given by other groups are in order of tens of keV, much smaller than the upper bounds set 
by recent experiments. The reason is obvious that the aforementioned reactions violate isospin 
conservation, there is a large suppression factor of about (1/43. 7) 2 , which reduces the widths 
by 3 orders. The calculations which are based on the assumption that the newly discovered 
-D*j (2317) and Z) s j (2460) are p-wave excited states of D s and D* predict their widths to be at 
order of a few to tens of keV. By contrary, if they are four-quark states, or molecular states, there 
may be more decay channels available, i.e. some modes are not constrained by the OZI rule, 
thus much larger total widths might be expected in that scenario. The authors of refs.[9l 12*0*] . 
for example, suggested that D*j(2317) and D s j(2460) are in non-cs structure (four-quark states 
etc.), and obtained much larger rates, even though still obviously smaller than the experimental 
upper bounds. So far, it is hard to conclude if they are p-wave chiral partners of D s and D* or 
four-quark states yet. 

We hope that the further more precise measurements of Babar, Belle and CLEO may offer 
more information by which we may determine the structure of the newly discovered mesons. 
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Appendix 



° ~ 2^ ' (19) 
h(-m v /2) - i 3 K/2) + u[I 3 (A s ) - I 3 (A H )] _ Q[A S - a;m„/2] 

1 2m v (l-uJ 2 ) 1-uj 2 ' 1 ' 
n -h(A s ) + I 3 (A H ) - m[I 3 (-m v /2) - I 3 (m v /2)] OK/2 - A s uj] 

° 2 ~ 2m v (l-u 2 ) T^? ' (21) 

B l 2loB 4 -B 2 -B 3 

° 3 = T + 2(1 -u;^ ' (22) 
n -Bi , 3B 2 -6u;B i + B 3 (2uj 2 + l) 

UA ~ 2(1 + 2(1-^)2 > W 



U5 ~ 2(1- W 2 ) + 2(1 -oj 2 ^ 2 



Bicu 2B 4 (2co 2 + l)-3co(B 2 + B 3 ) 

° 6 ~ 2(l-^) + 2(1-^)2 > ( 25 ) 

fix = m 2 0-J 3 (A H ), (26) 

Z3 2 = A s O 1 — '■ (um v + 2A s ), (27) 



*3 = f + ^W + 5 As/3(As) _ Aff/3(A4 



^ = m n A s O | A s [J 3 (A s )-J 3 (A g ) | J 3 (m„/2) - I 3 (-m v /2) ^ 



2 2m r] 



/5(a2 ' a3 ' a4) = Jo dX l + 2x2(1 -a 4 ) + 2x(a 4 - 1) { 16^7* L 2 ^-pH^"^)] 



x[l + erf(e v ^)] + I ^^ ^exp[-y(m2-2e 2 )]}, (30) 
a 2 (l - x) + a 3 x 



^1 + 2(a 4 - l)x + 2(1 - a 4 )x^ 



^ = 16^72 ^ ^i/r^ ^exp[ s A 2 (x)][l + erf(A(x)v^)], (32) 



(31) 



167T 3 / 2 J1/A2 

n2 

I I ' 

|q| 

n 2 = (34) 



Til = H (33) 
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= fh, (35) 
|q| 

n 6 = (36) 
|q| |q| 

n 3 -n 3 -Tl 6 2K 5 

l<-4 — j— T5 j— T) l-joj 

|q| |q| 

A(x) = A s -*|q|, '2 = ^(0,^,^) (39) 

m = -l 2 -A s n, (40) 

n 2 = ^[-I 3 (A S ) + I 3 (A H )}, (41) 

7i 3 = A|^+(M + A 5 )/ 2 , (42) 

n 4 = ^-[A S I 3 (A S ) - A H I 3 (A H )}, (43) 

n 5 = ^[I 3 (A S ) - I 3 (A H )}, (44) 

?i 6 = m 2 s K-I 3 (A H ), (45) 
where q is the three momentum of the emitted photon. 



References 

[1] The Babar Collaboration, B.Aubert et al., Phys. Rev. Lett. B 90, 242001 (2003); F. Porter, 
Eur. Phys. J. C33, S219-S222 (2004). 

[2] The CLEO Collaboration, D. Besson et al., Phys. Rev. D68, 032002 (2003). 

[3] Belle Collaboration, P. Krokovny et al., Phys. Rev. Lett 91, 262002 (2003). 

[4] W.A. Bardeen, E.J. Eichten and C.T. Hill, Phys. Rev. D 68 , 054024 (2003); M.A. Nowak, 
M. Rho and I. Zahed, Acta Phys.Polon. B35, 2377-2392 (2004) ;E. Kolomeitsev and M. 
Lutz, Phys.Lett.B582, 39 (2004). 

[5] K.D. Chao, Phys. Lett. B599, 43-47 (2004). 

[6] S. Narison, Phys. Lett. B605, 319-325 (2005). 

[7] E.V. Beveren, G. Rupp, Phys. Rev. Lett. 91,012003 (2003); E.V. Beveren, G. Rupp, Eur. 
Phys. J. C32, 493-499 (2004). 

[8] Y.Q. Chen and X.Q. Li, Phys. Rev. Lett. 93, 232001 (2004); A.P. Szczepaniak, Phys. Lett. 
B567, 23-26 (2003); E. Beveren and G. Rupp, |hep^ph/0305035| T.E. Browder, S. Pakvasa 
and A. A. Petrov, Phys. Lett. B578, 365 (2004). 



11 



[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 

[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 



H.-Y. Cheng and W.-S. Hou, Phys. Lett. B566, 193 (2003). 

T. Barnes, F. Close and H. Lipkin, Phys. Rev. D68, 054006 (2003). 

A.P. Szczepaniak, Phys. Lett. B567,23 (2003). 

Xin-Heng Guo, Hong-Wei Ke, Xue-Qian Li, Xiang Liu and Shu-Min Zhao, arXiv: 
hep-ph/0510146. 



M. Nielsen, arXiv: hep-ph/0510277 



Wei Wei, Peng-Zhi Huang and Shi-Lin Zhu, arXiv: hep-ph/05 10039 



P. Colangelo and F. De. Fazio, Phys. Lett. B570, 180 (2003). 

W.A. Bardeen, E.J. Eichten and C.T. Hill, Phys. Rev. D68, 054024 (2003). 

S. Godfrey, Phys. Lett. B568, 254 (2003). 

Fayyazuddin adn Riazuddin, Phys. Rev. D69, 114008 (2004). 

P. Colangelo, F. De Fazio and A. Ozpineci, Phys. Rev. D72, 074004 (2005). 

S. Ishida et al., arXiv: hep-ph/0310061 



The SELEX Collaboration, Phys.Rev.Lett. 93 242001 (2004). 



Babar Collaboration, arXiv: hep-ex/0408087 



Belle Collaboration, B. Yabsley, arXiv: hep-ex/0507028 



FOCUS Collaboration, R. Kutschke, E831-doc-701-v2, available at www-focus.fnal.gov. 
A.D. Polosa, Riv. Nuovo Cim. 23N11, 1, 2000. 

D. Ebert, T. Feldmann and H. Reinhardt, Phys. Lett. B388, 154-160 (1996); T. Feldmann, 
arXiv: [hep-ph/9606451; D. Ebert, T. Feldmann, R. Friedrich, H. Reinhardt, Nucl. Phys. 
B434 619-646 (1995). 



[27] A. Deandrea, N. Di Bartolomeo, R. Gatto, G. Nardulli and A.D. Polosa, Phys.Rev. D58, 
034004 (1998); A. Deandrea, R. Gatto, G. Nardulli and A.D. Polosa, Phys.Rev. D59, 
074012 (1999); A. Deandrea, R. Gatto, G. Nardulli and A.D. Polosa, J. High Energy Phys. 
02, (1999) 021; A. Deandrea, R. Gatto, G. Nardulli and A.D. Polosa, Phys.Rev. D61, 
017502 (2000); A.D. Polosa, arXiv: |hep-ex/9909371| 

[28] A. Deandrea, G. Nardulli and A.D. Polosa, Phys.Rev. D68, 097501 (2003). 

[29] S. Okubo, Phys. Lett. 5, 165 (1963). G. Zweig, CERN Rep. 8419/TH-412; CERN Preprints 
TH-401, TH-412. J. Iizuka, Prog. Theor. Phys. Suppl. 37/38, 21 (1966). 

[30] P. Cho and M.B. Wise, Phys. Rev. D49, 6228 (1994). 

[31] Belle Collaboration, Y. Mikami et al., Phys. Rev. Lett 92, 012002 (2004); 



12 



[32] Babar Collaboration, B. Aubert et al., arXiv: hep-ex/0408067 



[33] F. Giacosa, Th. Gutsche, V.E. Lyubovitski and Amand Faessler, Phys. Rev. D 72, 094006 
(2005). 

[34] The Data Group, Phys. Lett. B592, 1 (2004). 

[35] J. Gasser and H. Leutwyler, Nucl. Phys. B250, 465 (1985). 



13 



